High-energy behavior of hadronic total cross sections from lattice QCD 
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Abstract 

By means of a nonperturbative approach to soft high-energy hadron-hadron scattering, based on the analytic continua- 
tion of Wilson-loop correlation functions from Euclidean to Minkowskian theory, we shall investigate the asymptotic 
energy dependence of hadron-hadron total cross sections in lattice QCD: we will show, using best fits of the lattice 
data with proper functional forms satisfying unitarity and other physical constraints, how indications emerge in favor 
of a universal asymptotic high-energy behavior of the kind B log 2 s for hadronic total cross sections. 
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1. Introduction 

Among the oldest open problems of hadronic physics, 
not yet satisfactorily solved in QCD, there is the prob- 
lem of predicting hadronic total cross sections at high 
energy from first principles. Present-day experimen- 
tal observations (up to a center-of-mass total energy 
yfs = 7 TeV, reached at the LHC pp collider |l]]) seem to 
support the following asymptotic high-energy behavior: 
cr*'(s) ~ Slog 2 i, with a universal (i.e., not depending 
on the particular hadrons involved) coefficient B ^ 0.3 
mb J2]. This behavior is consistent with the well-known 
Froissart-Lukaszuk-Martin (FLM) theorem y|], accord- 
ing to which, for s — » oo, crf^{s) < (7r/m 2 ) log 2 (s/so), 
where m n is the pion mass and sq is an unspecified 
squared mass scale. As we believe QCD to be the fun- 
damental theory of strong interactions, we also expect 
that it correctly predicts from first principles the behav- 
ior of hadronic total cross sections. However, in spite 
of all the efforts, a satisfactory solution to this problem 
is still lacking. (For some theoretical supports to the 
universality of B, see Ref. [4] and references therein.) 

This problem is part of the more general problem 
of high-energy elastic scattering at low transferred mo- 
mentum, the so-called soft high-energy scattering. As 
soft high-energy processes possess two different energy 
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scales, the total center-of-mass energy squared s and the 
transferred momentum squared t, smaller than the typi- 
cal energy scale of strong interactions (|f| < 1 GeV 2 «: 
s), we cannot fully rely on perturbation theory (PT). 
A nonperturbative (NP) approach in the framework of 
QCD has been proposed in ^ and further developed 
in |0] : using a functional integral approach, high-energy 
hadron-hadron elastic scattering amplitudes are shown 
to be governed by the correlation function (CF) of cer- 
tain Wilson loops defined in Minkowski space [60- This 
CF can be reconstructed by analytic continuation from 
the CF of two Euclidean Wilson loops that can 

be calculated using the NP methods of Euclidean Field 
Theory. The analytic-continuation relations have al- 
lowed the NP investigation of CFs using some ana- 
lytical models, such as the Stochastic Vacuum Model 
(SVM) Jul, the Instanton Liquid Model (ILM) fTH[l2|l , 
the AdS/CFT correspondence 11311 . and they have also 
allowed a numerical study by Monte Carlo simulations 
in Lattice Gauge Theory (LGT) 1U21I14I1 . 

In what follows, after a brief survey of the NP ap- 
proach to soft high-energy scattering in the case of 
meson-meson elastic scattering, and of the numerical 
approach based on LGT, we will focus on the search for 
a new parameterization of the (Euclidean) CF that, in 
order: i) fits well the lattice data; z'z) satisfies unitarity 
after analytic continuation; and, most importantly, Hi) 
leads to a rising behavior of total cross sections at high 
energy as B log 2 s, in agreement with experimental data 
B15I1 . In our approach, the coefficient B turns out to be 
universal, i.e, the same for all hadronic scattering pro- 
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cesses, being related to the mass-scale /j which sets the 
large impact-parameter behavior of the CF. 

2. High-energy meson-meson elastic scattering am- 
plitude and Wilson-loop correlation functions 

In the soft high-energy regime, the elastic scattering am- 
plitude M(hh) of two mesons, of the same mass m for 
simplicity, can be reconstructed from the scattering am- 
plitude M(dd) of two dipoles of fixed transverse sizes 
f\,2x, and fixed longitudinal-momentum fractions of 
the quarks in the two dipoles, after folding with squared 
wave functions p lt2 = |^i,2l 2 describing the interacting 
hadrons 01 , 

M(hh){s,t) = J d 2 v pi(vi)p2(v 2 )M(dd)(s,C, v\, v 2 ) 

= {{M {d d){s,t-\,2))), (1) 

where v, = (fix,fd denotes collectively the dipole vari- 
ables, d 2 v - dv\dv2, J dvi = J d 2 fx JT dfi, and 
J dvj pi(vi) — 1. In turn, the dipole-dipole (dd) scat- 
tering amplitude is obtained from the (properly normal- 
ized) CF of two Wilson loops (WL) in the fundamental 
representation, defined in Minkowski spacetime, run- 
ning along the paths made up of the quark and anti- 
quark classical straight-line trajectories, and thus form- 
ing a hyperbolic angle y — log(s/m 2 ) in the longi- 
tudinal plane. The paths are cut at proper times +T 
as an infrared regularization, and closed by straight- 
line "links" in the transverse plane, in order to en- 
sure gauge invariance; eventually, T — > oo. It has 
been shown in |7H9[] that the relevant Minkowskian 
CF @m(x> T\ V\,V2) (Zx being the impact parameter, 
i.e., the transverse separation between the two dipoles) 
can be reconstructed, by means of analytic continua- 
tion, from the Euclidean CF of two Euclidean WL, 
&e(0;T-J ± ;v u v 2 ) = (W?'Wn)/(('W?)('W?)) - 1, 
where (. . .) is the average in the sense of the Eu- 
clidean QCD functional integral. The Euclidean WL 
f Wf 2 = N^TriTexpl-ig <f> Ci A^dXfJ] are calculated 
on the following quark [g]-antiquark [q] straight-line 
paths, Ct : Xf\r) = Zi + %t + ff'n, with r € [-T, T], 
and closed by straight-line paths in the transverse plane 
at r = +T. Here p\^ = m(± sin §, 0j_, cos |), r,- = 
(0,^,0), Zi = tffl(0,£,0) and f< = 1 - f* = -f t . 
We define also the CFs with the infrared cutoff removed 
as Ce,m = limbec Qe,m- The dd scattering amplitude is 
then obtained from Ce(9', ■ . .) [with 9 e (0, n)] by means 
of analytic continuation as (f = -\q±\ 2 ) 

M idd )(s, t; vi , v 2 ) = - i 2s f d 2 Zi_S^C M ix\ z±.;v x , v 2 ) 



Choosing p\ 2 invariant under rotations and under the 
exchange /j — > 1 — /; (see Refs. 0), Ce can be 
substituted in ([TJ with the following averaged CF: 
C™ (6; n , v 2 ) = f%f fi{C £ (^; v, , v 2 ) + 
C E (0;z±;vi,v 2 ) + C E (9;zxWi,v 2 ) + C E (0;z x ;vi,v 2 }, 
where f ix = |^J(cos fa, sin fa), V; = Qf ix \,fi) and v ; = 
{—fix, 1 - fi). Similarly, one defines the Minkowskian 
averaged CF, ■ As a consequence of the (Euclidean) 
crossing-symmetry relations II 1611 . Ce(x — 9', Z±', vi, v 2 )= 
C E (9;z±; v\,v 2 )-C e {9;zx, v u v 2 ), C| ve is automatically 
crossing-symmetric, i.e., C a ^ e {n -&,...)= C™ e (9; . . .). 



3. Wilson-loop correlation functions on the lattice 
and comparison with known analytical results 

In Refs. ITLZi [Till two of us performed a Monte Carlo 
calculation of Ce in quenched QCD at lattice spacing 
a(J3 — 6) =s 0.1 fm, on a 16 4 hypercubic lattice. We used 
loops of transverse size a at angles cot 9 = 0, ± 1 , +2. The 
longitudinal-momentum fractions were set to = | 
without loss of generality 11211 . We studied the con- 
figurations z± || II fix Czzz"), zx -L n± II 
("zyy") in the transverse plane, and the averaged quan- 
tity ("ave") defined above, for loops at transverse dis- 
tances d = \z.x\/a - 0, 1,2. 

Numerical simulations of LGT provide (within the 
errors) the true QCD expectation for Ce', approximate 
analytical calculations of Ce have then to be compared 
with the lattice data, in order to test the goodness of the 
approximations involved. Ce has been evaluated in the 
SVM, M| = § e -T*scote + insects _ 1 (^J, in PT , 

C ( ^^K v coi 2 9 fi EE in the ILM, C£ LM) = 



| UJ,[12|], and, using the AdS/CFT correspondence, for 



planar, strongly coupled N - 4 SYM at large \zx\, 



£>(AdS/CFT) _ e -j-jL + k 2 cot g + K 3 cos 
E 



Scots 



- 1 11311. The coef- 
ficients K, = Kj(zx',vi,v 2 ) are functions of Zx and of 
the dipole variables fx,fi- The comparison of the lat- 
tice data with these analytical calculations, performed 
in Ref. [14] by fitting the lattice data with the corre- 
sponding functional form, is not fully satisfactory, even 
though largely improved best fits have been obtained by 
combining the ILM and PT expressions into the expres- 
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Zxe iqA C E (0 -> -ix\Zx\ v u v 2 ) . 



(2) 



sion C% m 

dependence of total cross sections, the above analyti- 
cal models are absolutely unsatisfactory, as they do not 
lead to Froissart-like total cross sections at high energy, 
as experimental data seem to suggest. Infact, the SVM, 
PT, ILM and ILMp parameterizations lead to asymptot- 
ically constant cr[™, while the AdS/CFT result leads to 
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4. How a Froissart-like total cross section can be ob- 
tained 

We will now introduce, and partially justify, new param- 
eterizations of the CF that: i) fit well the data; ii) satisfy 
the unitarity condition after analytic continuation; and 
iii) lead to total cross sections rising as B log 2 s in the 
high-energy limit lll5ll . Regarding unitarity, from ([TJ 
and (f2]l one recognizes that the quantity A(s, \z±\) = 
((CmU 7 ; 2x\ 1, 2))) is the scattering amplitude in impact- 
parameter space, which must satisfy the unitarity con- 
straint \A + 1| < 1 (see ILLM]). Since j dvi p,(v,) = 1, 
this is the case if the following sufficient condition is 
satisfied (as we can replace Cm — > C a ^ when averaging 
over the dipole variables, a similar but weaker sufficient 
condition can be given in terms of C"^): 

|CAf(r,2x.vi,v 2 ) + l|<l vi,v 2 . (3) 

The conditions above constrain rather strongly the pos- 
sible parameterizations. For example, conditions ii) and 
iii) cannot be simultaneously satisfied if the angular de- 
pendence can be factorized, for in this case the unitarity 
constraint would imply c*'(y) — > const, for^- — > 00. 
We shall assume that the Euclidean CF can be written 
as Ce - sxpKE - 1, where Ke = Ke(6;z ± ;v\,V2) is 
a real function (since Ce is real [14]). This assump- 
tion is rather well justified: in the large-iV c expansion, 
C E ~ 0(1 /W 2 ), so that C E + 1 > is certainly satis- 
fied for large N c ; all the known analytical models sat- 
isfy it; the lattice data of Refs. Ill 2l 11411 confirm it. 
The Minkowskian CF is then obtained after analytic 
continuation: Cm - expAOw - 1, with Km(x', ■ ■ ■) = 
Ke(0 — > -ix\ • • •)■ At large x, Cm is expected to obey 
the unitarity condition Q, which in this case reduces to 
ReK M < Vzj.,vi,V2. 

For a confining theory like QCD, Ce is expected to 
decay exponentially as Ce ~ ffc^^ 1 ' at large \z±\, where 
fi is some mass-scale proportional to the mass of the 
lightest glueball (Mq ^1.5 GeV) or maybe to the in- 
verse of the so-called vacuum correlation length A vac 
(e.g., yu = 2/Avac in the SVM), which has been mea- 
sured on the lattice [20], both in quenched (A vac 0.22 
fm) and full QCD {A vac 0.30 fm). Therefore, we 
should require the same large-|zj behavior for Ke, i.e., 
K E ~ e"* 1 . 

Let us now assume that the leading term of the 
Minkowskian CF for \ ~ * +°° is of the form Cm ~ 
exp(/yS/Of)e _/J ' Zi ')-l [recall^ log(j/w 2 )], where B- 
P(v\ , V2) is a function of the dipole variables and f(x) is 
a real function such that fix) — > +°° for^ — » +00. In 
this case, the unitarity condition 0) is equivalent (for 
large x) to FmyS > 0. This z± dependence is expected 



to be valid only for large enough \z±\, but for simplic- 
ity we shall first assume that it is valid V|zJ > 0. By 
virtue of the optical theorem, cr'^(s) ~ s~ x lmM(hh)(s,t = 
0), we have that cr[™ ~ 47r/T 2 Re«./(77,j8)», where 
J{r,,B) = f °°dyy[l - exp(^ e «')], with fix) = e\ 
and y = /-i\z±\- Expanding the exponential, integrat- 
ing term by term, and deriving with respect to 77, we 
find dJ/drj = - Z7=i (-z)"/(n\n) = E t {z) + log(z) + y, 
for I arg(z)| < n (z — -iBe n ), where y is the Euler- 
Mascheroni constant and E\(z) is Schlomilch's expo- 
nential integral (see 12111 ). Since E\{z) ~ e~ z /z at 
large |z|, for Rez > o ImB > 0, the asymptotic 
form of dJ/drj is readily obtained; re-integrating in 77 
and substituting back 77 = log fix), we find ~ 
4^- 2 «i log 2 fix) + log /Cf)0og \B\ + y) + . . .». If one 
takes fix) = X Pe " x ^ me resulting asymptotic behavior 
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The same result is obtained assuming the above approx- 
imation for Cm only for \% ± \ > zo » pT x , |^j_|: the dif- 
ference in erf* , coming from the integration of Cm over 
the finite region < zo, is bounded by a constant due 
to the unitarity constraint. The analysis can be repeated 
for C ave without altering any conclusion. We want to 
emphasize that the above result is universal, depending 
only on the mass scale ft, which sets the large-|zxl de- 
pendence of the CF, since the integration over the dipole 
variables does not affect the leading term. The universal 
coefficient B is not affected by the masses of the scatter- 
ing particles: for mesons of masses wi_2, the rapidity 
becomes^ ~ log(^j-), which simply corresponds to a 
change of the energy scale implicitly contained in ©. 



5. New analysis of the lattice data 

We show now three parameterizations = exp K ( £ - 1, 
i = 1,2,3, that satisfy the criteria i)-iii) listed above, 
together with the corresponding estimate of the asymp- 
totic total cross section at high energy [15]. We focus 
our analysis on the averaged CF C ave , that is "closer" 
to the hadronic scattering matrix M{hh)- As C are is 
crossing-symmetric, so are our parameterizations. 

In order to parameterize Ke, a possible strategy is to 
combine known QCD results and variations thereof. We 
have then exponentiated the two-gluon exchange and 
the one-instanton contribution (i.e., the ILMp expres- 
sion), adding a term which could yield a rising cross 
section, e.g., a term proportional to cos 6 cot 6, as in the 
AdS/CFT result. We thus find the following parameter- 
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Table 1 : Chi-squared per degree of freedom for a best fit with the indicated 
function. 
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Figure 1 : Lattice data for £75* and best fit with Cjf . 

strategy is suggested again by the AdS/CFT result: one 
can try to adapt to QCD analytical expressions obtained 
in related models, such as N = 4 SYM. Although QCD, 
of course, is not N = 4 SYM, it is sensible to assume 
a similar dependence on 6 (basically assuming the exis- 
tence of the yet unknown gravity dual for QCD). In this 
spirit, the second parameterization that we propose is: 
K e = Jr7j +K 2(% ~ 0) cot 6 + K 3 cos cote. Beside the 
AdS/CFT-like terms, it contains also a 6 cot 9 term. Out- 
last parameterization is: = + Kj_{\ — &) 3 cos 8. 
While the first term is "familiar", the second one is not 
present in the known analytical models, but it is a fact 
that the resulting best fit is extremely good (see Fig.[TJ. 
In Table Q] we report the values of the chi-squared per 
degree of freedom (x\ f ) of the best fits to the lattice 
data. 

In the three cases, the unitarity condition RcK^ < 
is satisfied if K2 > 0: this is actually the case for our best 
fits (within the errors). The leading term after analytic 
continuation is of the iormx p e x which, according to 
leads to cr™ ~ filog 2 s. The value of B = 2n/n 2 , ob- 
tained through a fit of the coefficient of the leading term 
with an exponential function ~ <? over the available 
distances, is found to be compatible with the experimen- 
tal result (within the large errors) in all the three cases 
(see Table 0. However, this must be taken only as an 
estimate, as lattice data are available only for small IzJ. 

6. Conclusions 

We have shown how a universal and Froissart-like 
hadron-hadron total cross section at high energy can 
emerge in QCD, and we have found indications for 
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Table 2: Mass-scale fi, "decay length" A = l/fi and the coefficient B — 2tt//j 2 
obtained with our parameterizations. 

this behavior from the lattice. The functional integral 
approach provides the "natural" setting for achieving 
this result, since it encodes the energy dependence of 
hadronic scattering amplitudes in a single elementary 
object, i.e., the loop-loop CF. 
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